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Fluctuation and dissipation dynamics is examined at all temperature ranges for the general case 
of a background time evolving scalar field coupled to heavy intermediate quantum fields which in 
turn are coupled to light quantum fields. The evolution of the background field induces particle 
production from the light fields through the action of the intermediate catalyzing heavy fields. Such 
field configurations are generically present in most particle physics models, including Grand Unified 
and Supersymmetry theories, with application of this mechanism possible in inflation, heavy ion 
collision and phase transition dynamics. The effective evolution equation for the background field is 
obtained and a fiuctuation-dissipation theorem is derived for this system. The effective evolution in 
general is nonlocal in time. Appropriate conditions are found for when these time nonlocal effects 
can be approximated by local terms. Here careful distinction is made between a local expansion 
and the special case of a derivative expansion to all orders, which requires analytic behavior of the 
^ ^' evolution equation in Fourier space. 
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p ,■ I. INTRODUCTION 

CIh' . . . . 

^ ' Dissipation physics can play an important role in systems governed by quantum fields in subject areas as diverse 
, as condensed matter to particle physics and early universe cosmology. Dissipation effects can have experimental 
consequences in many systems. In cosmology, inflation could occur concurrently with particle production due to 
' dissipation from the scalar inflaton field in the context of the warm inflation scenario fr, leading to observable 
^ signatures in density perturbations [3, In heavy ion collisions, the chiral symmetry QCD transition could be 
Cf^ ' accompanied by dissipative effects Q. More generally, any phase transition in a condensed matter, particle physics or 
0^ cosmological system will be accompanied by dissipative effects during the dynamical evolution of the order parameter 
C characterizing the phase transition ^61] . 

The flrst principles origin of dissipation from quantum fleld theory has been a subject of interest for several decades. 
\^ The typical system studied 0, H, H, [l^j [HI fisl. Il4| is where a background scalar fleld is coupled to some variety of 
' bosonic and fermionic light quantum fields, with the entire system starting off in some statistical state. The general 
outcome is the time evolution of the scalar fleld induces production of particles from the light quantum flelds, creating 
dissipative and fluctuating effects. Such dynamics is of interest in particle physics, since this direct dissipative 
mechanism arises from very generic coupling schemes and also because such dynamical situations commonly present 
themselves in problems such as inflationary cosmology, heavy ion physics and in phase transitions in general. 
^ I , In particle physics models there is a novel interaction conflguration which also has been shown to lead to considerable 
' dissipative effects, and that is the case where a background scalar field is coupled to heavy intermediate quantum 
fields which in turn are coupled to the light quantum flelds [l^, [l^ . In this two-stage mechanism, the background 
scalar fleld indirectly induces particle production in the light flelds through the intermediate heavy fields which in a 
sense help to catalyze the effect. Dissipation from such field configurations in systems near thermal equilibrium has 
been shown to arise [l3| and application of these effects to inflation have been shown to have signiflcant importance 

A full understanding of this mechanism requires understanding not only the dissipative effects but also the fluctu- 
ations created and the relation between the two in the context of a fluctuation-dissipation theorem. This so far has 
not been done for this two-stage mechanism and will be one of the results of this paper. 

To arrive at this result, what will be required is to obtain the effective evolution equation for the background scalar 
field once all other fields to which it is coupled are integrated out. The resulting effective evolution equation, which 
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will be derived using the Schwinger closed-time formalism [19|, in general will have time nonlocal terms. The second 
interest of this paper is in determining when such nonlocal terms can be well approximated by local terms. In general 
this requires a careful examination of the different time scales involved in the dynamics of the system in study to 
determine the importance of nonlocal effects. For instance, when the characteristic time scale of the nonlocal (or 
memory) terms in the effective equation of motion is close to the typical reaction time of the system, it has been 
shown that large effects on the dynamics can arise as a consequence of the nonlocal terms as compared to the results 
coming from a local approximation for the equations [29|. In this paper we adopted the general approach of directly 
analyzing the nonlocal terms and then determining the relevant time scale for when we can approximate them as local 
forms. 

In practical application of these types of equations, there is an immense saving of effort as well as much more 
transparent understanding of the physics from a local equation as opposed to a nonlocal one, since the former can 
generally be analyzed with much less numerical treatment than the latter. Thus it is a very important question when 
and how accurately the generally nonlocal effective equations can be approximated by a local form. One very special 
type of local approximation is the derivative expansion 0, [HI [H, 0, El, US]- This is just a Taylor expansion of the 
nonlocal terms about some time point. Such an expansion therefore requires analytic behavior of the relevant nonlocal 
terms. However that may not arise in all circumstances, and for example even simple logarithmic corrections would 
make the derivative expansion invalid. Nevertheless, a careful distinction must be made between localized behavior 
as opposed to localized analytic behavior. Even in the general circumstance where the behavior is not analytic, the 
nonlocal terms may still be fairly peaked functions that can be well approximated as local. The nonlocal terms 
in the effective evolution equation we derive for this two-stage mechanism will be examined for when they can be 
approximated as localized. Although in this paper we carry out this investigation in the context of the two-stage 
mechanism, most of our significant results here would also apply to the direct dissipative mechanism (for a thorough 
discussion about the two-stage decay mechanism contrasted to the direct one, see e.g. Ref. [11,13). 

We assume throughout that the system is enclosed by an environment which is initially close to thermal equilibrium 
and remains in that state due to the existence of a relatively short interaction timescale in the various components 
contributing to the environment. In this respect we have a different situation than the one considered in the original 
work on the two stage mechanism [ill . [l6j and more recent numerical work [23]. In fact, we find that the local 
dissipative effects cease in the limit that the temperature T — > 0, as one might expect and in agreement with [itI. |23|. 

This paper is organized as follows. In Sec. II we give the Lagrangian density model used to study the two-stage 
decay mechanism and the respective effective equation of motion for a background scalar field that exhibits both 
the characteristic dissipation and fluctuation (stochastic) like terms. In Sec. Ill we derive the various generalized 
fluctuation-dissipation relations for the model under study and show the generality of these relations. In Sec. IV we 
define the local limit for the noise and dissipation terms and then we provide a thorough study both analytically and 
numerical to support localization in time for the noise and dissipation kernels appearing in the effective equation of 
motion given in Sec. II. Finally, our concluding remarks are given in Sec. V. 



II. THE MODEL 



We consider the following model of a scalar field <I> interacting indirectly through a catalyst field x to a scalar field 
a, which we can think of as our heat bath. The field $ has no direct coupling terms to the field a. The Lagrangian 
density for this model is the following 



1 TTl \ 

1 JJ-j 2 

+ -d^ad^a ~ - hMxcT^ (2.1) 

where is a parameter with mass dimension added in order to make the coupling constant h dimensionless. 

We consider $(a;) = ipdx) + (f>{x) where ipc{x) is a background like component of the scalar field, or vacuum 
expectation value of <f>, and (l){x) are the quantum fluctuations around this background field. The nonequilibrium 
equation of motion for ipc can be derived using for example the functional closed time path (CTP) formulation of 
Schwinger-Keldysh [l9j . In the CTP formalism the time integration is along a contour c from — cxd to +qo and then 
back to — oo. For reviews please see, for example, Refs. [13, IH, HI] . The detailed derivation of the effective equation 
of motion (EOM) for the field (pc was explained in many previous papers (see e.g. [HI, HB, El! ISS] ) ^^'^ therefore we 
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will only give here the final result for it. The final result is expressed in terms of effective propagators (that already 
include the effects of the bath fields). In the CTP formalism these field propagators, e.g. for (f), are given by 
(which also apply to any of the other boson fields): 

G++{x,x')=i{T+c^{x)c^{x')) , 
G^^{x,x) = i{T^(j){x)(j){x')) , 
G+-{x,x')^tWx')^{x)) , 

G^+(x,x')=*(0(x)0(x')), (2.2) 

where T+ and T_ indicate chronological and anti-chronological ordering, respectively. Gj"'' is the usual Feynman 
propagator. The other three propagators come as a consequence of the time contour and are considered as auxiliary 
(unphysical) propagators. From (|2.2p it also follows that the boson field propagators can also be expressed as 



G++{x,x') = G>{x,x')9{t - t') + G<{x,x')e{t' - t), 
G^- {x, x') = G> {x, x')e{t' -t)+G< {x, x')e{t - t'), 
G+~{x,x') = G<{x,x'), 
G^+{x,x')^G>{x,x'). 



(2.3) 



In deriving the effective EOM for ipc we first integrate over the a field. In this first step the Green's functions 
for the fields (in this case for the x field in leading order) coupled to the fermions get dressed by a field corrections. 
Subsequently functional integration over the other fields (in our case the x </> field fluctuations) is performed. In 
this process, we can express the boson Green's functions entering in the ipc effective EOM in the generic form 



G{x,x')+ / ct^zT,{x, z) G{z, x') — 5{x, x') , 



(2.4) 



where G{x, x') and S(x, x') are the (matrix) propagator and self-energy, respectively. The mass term in the above 
equation is the one appropriate the the respective field. 

In terms of a loop expansion for the fluctuation field (f) and x ^^'^d the generic expressions for the field propagators 
given above, the effective equation of motion that emerges for is of the form J^, 



+V,{x) / d^x'ifKx') 



-Im 



Gl+{x,x') 



25^Im 



G++(x,x') \e{t-t') 



d x'(pc{x') 



A2 



-Im 



G++{x,x') 



AgHm 



G++ (x, x')G++ [x, x')G++ (x, x') \0{t- t') 



Vc{x)£,i{x) +£,2{x) 



(2.5) 



where denotes the renormalized effective mass for </i that includes all local contributions. The nonlocal contribu- 
tions can be associated to dissipation and noise terms. For instance, the third and forth terms in Eq. (|2.5p are related 
to dissipation terms. The source terms and ^2 are stochastic (noise) fields associated with the imaginary terms in 
the effective action coming from the real-time evaluation of the perturbative loop diagrams leading to Eq (|2.5p . They 
have two-point correlation functions given by 



(a(x)ei(a:')) = yRe 



G++{x,x') +2/Re[Gj+(x,x')]' , 



(2.6) 



and 



{Ux)i2{x')) ^ -Re 



Gl+{x,x') 



2g^Re 



G++{x,x'fGl+{x,x') 



(2.7) 
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and are colored (space-time dependent) in general and Gaussian distributed, with probability distributions given by 
{Ni and are appropriate normalization constants) 



P[a] = iVf' exp . 



-Re 



G++{x,x') 



3 



Gt^{x,x') 



6(^') 



(2.8) 



and 



P[^2\ = ^2^' exp I -± / d^xd^x'Ux) 



-Re 



G++(x,x') 



+ 254Re 



G++{x,x'fG++{x,x') 



U^') • (2.9) 



Both dissipation and noise kernel are related to generalized fluctuation-dissipation relations, as we show next. 

III. FLUCTUATION-DISSIPATION RELATIONS 

A generalized fluctuation-dissipation relation between the nonlocal kernels in Eq. (|2.5p and the stochastic flelds 
correlation functions were derived in Ref. [11], while explicit expressions derived from the dissipation and noise kernels 
appearing in (|2.5p were obtained by Yokoyama in Ref. [131 ■ General expressions that relate both kernels can also 
easily be derived in the following way. Lets first define dissipative kernels, I?,;, as given by 



where the Ci=i,2 are defined by 



r n 2 r n 2 

Ci(x-x',i-i') = A^Im G++{x,x') sgn(t - t') H- 4/lm G++{x,x') sgn{t-t'), 



(3.1) 



(3.2) 



and 



C2(x - x',t - t') 



-Im 



G++{x,x') sgn{t - t') + 4gHm G++{x,x')G++{x,x')G++{x,x') sgn(i - t') , (3.3) 



where sgn{t — t') = 9{t — t') — 9{t' — t) is the sign function (note that only the part for t > t' oi these functions enters 
in the equation of motion Eq. (|2.5p ). The extra factor 2 in the definition of Eq. (|3.2p is for convention purposes, such 
that both Ci=i,2 will satisfy the same fluctuation-dissipation relation. Using Eqs. (|3.2p and p.3p in (|2.5p it follows 
that 



[U + Ml] ^,[x) + -^l{x) 



ipc{x) 



J d^x'ipl{x',t)Vi{x-x',0) ~ J d^xVc(x',<)P2(x-x',0) 



+ipc{x) I d^x' I dtVc(x',t')<y5c(x',t')2'i(x-x',i-t') + / d^x' f diVc(x', t')2'2(x - x', t - i') 

= ^c{x)ii{x) + i2{x) . (3.4) 



The third and forth terms in Eq. ()3.4|) can be shown to correspond, in the homogeneous approximation for the 
fleld to the one- loop $ vertex correction and the two- loop setting sun like mass corrections, respectively, for the 
model (I2.ip . Its is clear now from p.4p that Vi are dissipative like terms. That they satisfy generalized fluctuation- 
dissipation relations with the noise kernels defined by Eqs. p.6p and (|2.7p can be demonstrated when space-time 
Fourier transforming both kernels. We can write the Fourier transformed kernel JF[Ci(x — x', t — t')] = Ci{p, uj), where 
the Fourier transform J-'[- ■ ■] is defined by 
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^[■■■] ^ j Sxdt\ 



-ip-x+icji 



(3.5) 



Note that, from the Eqs. (|3.2p and (|3.3p . the kernels Ci are anti-symmetrical in the time. It then follows that their 
time Fourier transforms are purely imaginary. Thus, we can write \27\ 



where ri(p,a;) is a real quantity. The transformed dissipation kernels become 



(3.6) 



p,(p,w) = 2r,(p,c^) 



(3.7) 



Similarly, for the noise terms we introduce the kernels M(x — x', i — t'), where from Eqs. (|2.6p and (|2.7p . they are 
defined by 



and 



AAi(x-x',<-i') = YRe[G++(a;,a;')J + 2/Rc [Gj+(x, x')] , 



(3.8) 



A/'2(x-x',t-i') = ^Re [g++{x,x')^ 



G^+[x,x'fGl+[x,x') 



(3.9) 



and using again the properties of the field propagators in the CTP formalism, we notice that Mi are symmetrical in 
both space and time. Thus, their Fourier transforms in space and time, J-[Afi{x, x')] = Afi{p, w), are purely real. 

From this point on, we will restrict ourselves to parameter values such that A — 0{g'^) and ^ 1. In this situation 
we can drop the first term in the right-hand-side of both Eqs. (|3.8p and p.9p and likewise in the expressions for the 
kernels p.2p and (|3.3p . Repeating the analysis with these terms included presents no specific difficulties, but it is not 
particularly illuminating. 

In order to relate the dissipation and noise kernels, it is convenient to write the Feynman propagator G^^ (x, x') in 
real and imaginary parts as 



G++ {x, x')=F{x,x') + -p{x, x')sgn{t - t') 



(3.10) 



where F{x, x') and p{x, x') are the anti-commutator and commutator (or the statistical and spectral functions) of two 
fields defined, respectively, by 



F{x,x') = '-mx)A{x')}) = \ [G<{x,x') + G>{x,x')\ 



(3.11) 



and 



p{x,x') = imx),(^{x')]) = [G>{x,x')-G<{x,x')\ 



(3.12) 



For real fields, F{x,x') and p{x,x') are real and pure imaginary functions respectively, and have the symmetry 
properties F{x,x') = F{x',x) and p{x,x') = —p{x',x), which easily follow from the definitions (|2.2p and (|2.3p . Note 
also that since p{x,x') is an anti-symmetrical function, its Fourier transform, !F[p{x, x')] = p, is real. 
In terms of (|3.1ip and (|3.12p , the noise kernel A/i becomes 



A/'i(x-x',i-t') = 2g'^Re[G++{x,x')]' 



29' 



F^{ic - x', i - t' f + -p^i^ - x', i - t'f 



(3.13) 
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or in terms of Fourier transforms, 



^[A/-i(x-x',i-t')] = M(p, 



= 



UJ 



d^k f°° duj' 



Fx(k, uj')F^ (p - k, w - uj') + (k, c^')/5x(P - k, w - w') 



(.3.14) 



If we now use the property valid at equilibrium, coming from the Kubo-Martin-Schwinger relation for the propa- 
gators, G'^(p,w) = exp(/3cL')G'^ (p, w), it follows, from the definition of the spectral density, that [1^ 



G>(p,c^) = [l + n(w)]p(p,^) 
G<(p,cj) = n(tj)/5(p,w) , 



where n(u;) is the Bose-Einstein distribution function. 



•n{uj) 



1 



_ 1 ' 

and then we obtain for the Fourier transform of the statistical function, J^[F{x,x')] = F(p,aj), the result 

F(p,w) = i[l + 2n(L^)]p(p,u;) . 



Using Eq. (|3.17p and the identity 

1 + [1 + 2n{uj')][l + 2n(w - w')] = 2 (e'^" + l) n{Lj')n{Lj - uj') , 
we obtain for Eq. p.l4p the result 

/d^k dt,]' 
Wf y ^^"(^>('^-^')/5x(k,a;')/5x(p-k,c.-c.') 

Analogously, for the kernel Ci(x — x', f — i'), we find 

Ci(x-x',i-<') = 4/lm [Gj+(x,a;')]%gn(t-<') 

and expressing it in terms of transforms, after some algebra, we find 

/d^h r°° dt / 
(2^ j ^'^('^'W^-^')Px(k,^')Px(p-k,c.-c.') 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



Using now Eq. (|3.6p we finally find that the dissipation and noise kernel satisfy the generalized fluctuation-dissipation 
relation 



A/'i(p, cj) = 2tj n(ti;) + i 

Likewise, for the noise kernel M2 we flnd 

AA2(x-x',t-^') = 2g^Re\G++{x,x'fG++{x,x') 
= 2/ 



ri(p,c^) 



(3.22) 



1 



F^{^ - x', t - t'f + -p^i^ - x', f - t'f 



F^{^-^,t-t') 



iF^(x - x', i - <')Px(x - x', i - <')p^(x - x', < - 



(3.23) 
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and in terms of Fourier transforms, 

K~r , N 4 / fl, , .\ f d^k <Pq f°° duji du!2 / N / x / N 

X p.^{'k,uji)p^{q,uj2)pct,{p-k-ci,uj-uji-uj2) , (3.24) 
while for the transform of the kernel C2(x — x', < — t') we find 

C.(p,c.) = -2^g {e^ - l) y __„(..,)nMn(a; - a;i - c,) 

X p^(k, wi)p^(q, W2)p0(p - k - q, w - - W2) • (3.25) 

Upon using Eq. (|3.6p . we again find that the dissipation and noise kernels r2(p,a;) and 7V2(p,(x') are also related in 
an identical way to Eq. ()3.22p . 

These fluctuation-dissipation relations are actually generic to any product of two-point Green functions in the CTP 
formalism. To see this, begin by expressing the self-energy matrix elements in a similar way to the CTP propagators 
as 



J:—{x,x') = Y,> {x, x')9{t' - t) + J:< {x, x')9{t ~ t'), 
^+-{x,x') = -E<(x,a;'), 

^-+{x,x') ^ -T,>{x,x') , (3.26) 

From ()3.26|) the following property follows, 

S++(a;,x') -I- Y.+^{x,x') + Y.-+{x,x') + Y.— {x,x') = . (3.27) 

In addition, the self-energies in ()3.26p are found in general to satisfy conditions similar to those for the propagators, 
such as 



E>{x,x') = i:<{x',x) , 

= iS<(»(x,x') . (3.28) 

The real and imaginary parts for the self-energy component E++(p) in the CTP formalism can be separated 
by introducing the commutator and anticommutator functions, analogous to the ones used to express the causal 
propagator Eq. (|3.10p . 



iS>(<)(x,x') 



"(a;, x') = J:f{x, x') + ^^p{x, x') sgn{t - t') , (3.29) 



where 



Epix.x') ^ [^>{x,x') -Y.<{x,x')] , 

Ep{x,x') = ^[J:>{x,x') + J:<{x,x')] , (3.30) 

which are pure imaginary and real respectively. Turning to our earlier example, Eq. (|2.5p . we find that the dissipation 
kernels are related to Sp and the noise kernels to Sf- 

For the dressed propagators as defined in (|2.4p . we can express the transform of the spectral function, p(p,Ci;), in 
the form 

p(p,c^)=zL2-a;(p)2 + iEf + izSp) - t ( cu^ - u;{p f + 1-J:f - ^-itp] , (3.31) 
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where a; (p)^ — |pp+M^. The Breit-Wigner form of Eq. (I3.3ip suggests that we can introduce a generahzed relaxation 
time r(p, lu) of the particle state via the relation 

Ep(p,c^) =4co(p)t(p,c^)-i. (3.32) 
The equilibrium self-energies can be expressed in terms of Ep, 

E>(p,L^) - [l + n{u;)]Ep{p,Lu) , 

E<(p,cj) = n(cj)Ep(p,w) . (3.33) 
Using Eq. p.33p . we find a local fluctuation dissipation relationship between the Fourier transforms of S^? and Ep, 

(p, Lu) = ^[l + 2n{Lu)] Ep(p, Lu) , (3.34) 
where we have defined !F['E,p{x,x')] = E_F(p,Ct') and !F[Y.p{x,x')] — Ep(p,w). 

IV. LOCAL LIMIT FOR THE NOISE KERNELS 

We now study the possible localization of the effective EOM for the background field. We first note that nonlinear 
equations of the form of Eq. (|3.4p are impossible to solve analytically. They are also notoriously difficult to solve 
numerically, since they typically exhibit highly oscillating nonlocal kernels that can lead to errors which quickly build 
up and that are too hard to control, thus preventing any simple numerical attempt to solve these kind of equations. 
There are however different approaches that can be applied to analyze these equations. For instance, in those cases 
where the kernels exhibit a strong exponential damping in time, like in the analysis performed in Refs. [Tsl . [TB|, a 
numerical approach to analyze the nonlocal equation of motion for the background field is possible and these analysis 
have shown that a local approximation for the kernels, and then for the equation of motion, is in very good agreement 
with the full numerical solution of the nonlocal equation. There are however cases where such a strong damping 
behavior for the kernels are not possible (see e.g. Ref. ^T7|) and we must resort to other alternative analysis to 
determine how good a local approximation is for the nonlocal equation. For example what we are going to show 
below is the nonlocal kernel can be analyzed and used to infer the validity of the local approximation, in particular 
determine a time scale for locality; all this can be accomplished without reverting to solving complicated equations 
of motion like Eq. (|3.4p . This is our objective here. Note that since both dissipation and noise kernels are related, if 
we find that the noise kernels can be localized, then the same must apply to the dissipation kernels in the effective 
EOM. We shall consider homogeneous background fields, for which the main issue will be the localization in time. 
Throughout this section we assume that the relaxation time of the environment is sufficiently short to be neglected. 

First of all, we shall define what we mean by a localization of the noise kernels TV. The concept of localization 
depends on a particular context. For example, if we are interested in the origin of density fluctuations in the early 
universe we might be interested in localization of the noise kernel on timescales smaller than the Hubble time. We 
shall therefore introduce the concept of localization with respect to a particular timescale r. 

We shall say that a function f{t) is slowly varying on a timescale r if the Fourier transform /(lu) satisfies 

/(w) =0 for uj>2t:/t . (4.f) 

A kernel function IC(t) will be described as localized on a timescale r with accuracy f — e if the Fourier transform 
IC{uj) satisfies 



iCicu) - IC{0) 



< e for uj < lixjT . (4.2) 



/C(0) 

It follows that, for /(t) slowly varying on a timescale r and /C(t) localized on a timescale r, 



/oo 
/C(t - t')/(i')'^i' - ^(0)/(i) + 0(e) , 
-OO 



(4.3) 



9 



which is a formal way of stating that JC{t) « JC{0)5{t — t'). Note that, if ]C{lo) is analytic in w, then the kernel can 
be localized and the kernel admits a local derivative expansion. In general, the derivative expansion might not exist, 
even when the kernel is localized in the sense given above. 

To help understand the motivation for the above definition, it may be useful to think of a Gaussian function 
exp(— t^/r^^j) with width Ti/2- The corresponding width of the Fourier transform is 27r/Tx/2- By the definition 
adopted above, on the scale T1/2, we would say that the function is not fully localized but roughly speaking about 
half localized, thus localized with accuracy around 50%. For the Gaussian to appear localized in t-space, one actually 
should look at it from a much bigger scale of say lOri/2, which means a much smaller region in w-space near the 
maximum part of the gaussian where it does not vary too much, and thus e in Eq. (j4.2p would be small. 

As shown in the previous section, the kernels Cj are related to the dissipation kernels Vj by Eq. (|3.H) . which 
transformed using ui (we omit the spatial variable for homogeneous fields) become, 

P,(c.) = -C,H-2f,H. (4.4) 
However, from the fluctuation-dissipation relation (|3.22p , we have that 



2lj 



r,(^) 



(4.5) 



If a local approximation for the noise two-point correlation function exists, i.e. Af{t — t') — Nf)5{t — t'), then the noise 
amplitude A^o will satisfy 



No = 7^(0) = 2Tf(0) , (4.6) 

where we have used Eq. (14. 5p and taken the limit as w ^ 0. Note in particular from Eq. ()4.6|1 that there is no T = 
local limit for the noise correlation function and, thus, no local limit for the dissipation terms in the effective EOM 
(EH) at T = 0. 

Equation (|4.6p is one of our main results. This simple relationship between the noise and dissipation has appeared 
previously as a high temperature limit in various applications, but our derivation is quite general and depends only 
on the localizability of the noise and dissipation. We shall describe in the next subsection how, and under what 
circumstances, the localization in time is realized. 

The noise amplitude Af{0) can be determined as follows. The inverse transform of the one-loop and two-loop noise 
kernels Eqs. ([3?T9l) and ([3721]) is defined by 



.F-MM(p,c.)]=A/-,(x-x',i-i')- / 7^ r"-^expHP-(x-x') + zc.(i-<')]M(P,^)- (4.7) 

J (2TTy J^^ (Itt) 

For homogeneous fields, ifidx) = ^c{t), all loop diagrams in the effective action carries zero external space momentum, 
p = 0, and then immediately follows from Eq. (|4.7p that 7Vi(x — x', t — t') and A/2(x — x', t — t') are delta correlated 
in space (7\/i(x — x',< — t') ^ (5(x — x')). If also a local limit for the noise two-point correlations are valid, then the 
noise amplitudes become, from Eqs. (|3.19p and p.24p . respectively 

(2^7 ^ ^'^('^i)[l + "K)]/5x(k,^i), (4.8) 

while for the two-loop noise kernel. 



rr ^ d f d^k d^q dcui dLU2 , n / xr / m 

A/-.(o) = 2,^y (^(4y_^^"(-iw-^)[i+"(-i+-^)] 

X Px(k,wi)pj^(q,W2)p0(k-Kq,wi +W2) . (4.9) 

From now on we will assume the simplest case of homogeneous fields and study the localization for the noise kernels in 
time. If such a local approximation exists for the noise kernels, then Eq. (|4.6p suggests that the Fourier (w) transform, 
J\fj{uj), should not differ much from the lo = value, at least for some region in lu space, — ti^iocai ^ ^ ^5 '^locai. We 
now try to find this region in lo space. 
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A. The one-loop noise kernel 



For the general case, which should be valid for both the low as well high temperature limits, we can use the full 
expression for the spectral density, Eq. (|3.3ip . which in terms of the decay width = — ImSi/[4wi(k)], can be 
written in the form (i = 0, x) 



Pi(k,w) 



4w,(k)Tri(k,L^) 



[c.2-c.f(k)]^+[2^,(k)r-i(k,c.)]- 



(4.10) 



We will now study the localization of the one-loop noise kernel Eq. (|3.8p . Recall that here we are restricting our 
analysis to the parameter regime where A <^ g^, where the first term in Eq. (|3.8p can be neglected compared to the 



second one. The second term will only involve the spectral density for the x held. The decay width t~ receives 
contributions from both one and two a loop diagrams. Both results are well known (see e.g. Ref. [l3| for the first 
one and Ref. [l^ for the second). The single a loop diagram is 0{h?), while the two-loop diagram is 0{g'^,h'^). 
Therefore, for parameter values such that h ^ (these parameter values were in fact the ones shown to be relevant 
for the dynamics studied in Refs. [H, [3, EH [3 ) i the dominant contribution comes from the single loop. Thus, for 
the dominant decay channel we have that (fc = |k|) [l3| 



-i(k,c.) 



167rti;^(k) 



k)-e{- 



k)] 



In 



-f3\u)+k\/2 



-p\ui-k\/2 



which is valid for very light a bath fields, m„ « 0. 
Equation p.lOp then becomes (for p= |p| =0) 



(4.11) 



(Pk 
(2^ 



27 



n{uj')n{uj — w')p;^(k, a;')p^(k, w — to') 



(4.12) 



1. The low temperature limit 



At low temperatures T <C m^, we can easily verify that the dominant contributions to the noise kernel Eq. (|3.19[) 
come from ^ m^. We can approximate the spectral function Eq. (|3.3ip for the x field as 



4 



(4.13) 



with LOx^) ^ "'"X l|4.1ip . Thus, from Eqs. 

considered before, h f,g^ and X g^, we obtain 



p.l9p and (|4.13p . and for the same parameter values 



16g^ 



d^k 

(2^ 



27r 



(4.14) 



The dimensionful parameters can all be rescaled out of the integral to leave 



(4.15) 



In principle F only depends on /3a;, but in practice there is a breakdown of (|4.13p at w = fc when the light field mass 
Too- = 0, resulting in the h dependence. 

The results of a numerical evaluation of the function F(/i, /Jw) have been plotted in Fig. [U which takes the 
temperature T/m^ = 0.01 for the cases h = 0.1 (solid) and 1.0 (dashed). Below this temperature, the numerical 
calculations find almost no change in Y{h^[3u)), which is consistent with our parameterization. The best 3-parameter 
fit to the numerical curve at small values of a; = w/T is of the form 
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FIG. 1: The rescaled noise amplitude F(^,a::), where x = uj/T, at T/m^ = 0.01 for h = 0.1 (solid) and 1.0 (dashed). Below 
this temperature these curves do not change significantly. 



¥{h, x) = A{h) + B{h) X \n{x) + C{h) X , (4.16) 

where we find ^(0.1) = 0.0458, B{0.1) = 0.144 and C(O.l) = 0.126 and yl(l.O) = 0.0161, B(l.O) = 0.0068 and 
C(l.O) = —0.0004. These numerical results provide convincing evidence that the noise kernel can be localized in the 
sense described earlier. For example, we can say for h = 0.1 that the noise kernel is localized on a timescale 1600T~^ 
with 95% accuracy and for h = 1.0 the noise kernel is localized on a timescale 100T~^ with 95% accuracy. Note that, 
although this may be a large timescale compared to the scale of particle interactions, it may still be a small timescale 
compared to the timescale of a cosmological application. 

It is interesting that the numerical evidence supports the conclusion that the noise kernel is not analytic in lo. The 
logarithmic terms in u are related to subleading non-local effects. This also has an important implication for the 
dissipation terms. The dissipation kernels and 2?^ are related to the noise kernel by the fluctuation-dissipation 
relation Eq. (j4.5p . and in the low temperature limit, 

= -iujf)^ = -2i(3u}JV. (4.17) 

Therefore the dissipation terms can have a local analytic expansion at leading order to give a tp term, but beyond 
that there is not a local analytic expansion but rather terms of the form 

For a definition of the logarithmic derivative, see Ref. [30l |. 

2. The high temperature limit 

For high temperature T 3> m-^ and T 3> w, the specral function is concentrated in the range of values of lo' close 
to ±uj{p), where 

a;x(k) = (fc'+m2)i/2. (4.I9) 
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FIG. 2: The rescaled noise amplitude F{h,x), where x = uj/T, at T/rrix ~ 100.0 for h — 0.1 (solid) and 1.0 (dashed, increase 
by a factor 10,000). 



We take the distribution functions out of the energy integral. 



A/-i(u;)«-.9V- + l) j ^„(c.^(k))[l + n(c.^(k))] ^ ^ — p^(k, c.')/5;,(k, c - c') 
The spectral function is given by Eq. (|4.10p . with T^(fc) = T^(^,bj^{ky). After the w' integration, 



(4.20) 



(Pk AT^{k) n(tJx(k)) [1 + n{uo^{\^))] 



(27r)3 4 + Tx(fc)2tj2 



(4.21) 



Note that this result for the noise kernel is consistent with the earlier relationship A/i(0) = 2Tf (0), where f (0) is the 
high temperature dissipation coefficient first calculated by Hosoya and Sakagami [3] . 

The noise kernel is localized in the sense which we defined earlier. For a rough estimate, we introduce the relaxation 
timescale at k = T and take out of the k integral. The noise kernel is then localized at the timescale ~ 30r^ 
with accuracy 99%. Alternatively, we can say in less formal terms that the noise kernel appears localized for physics 
on timescales larger than the relaxation time. A numerical calculation of the noise amplitude is shown in Fig. [2] and 
it is consistent with these analytic estimates of the localization scale. 



B. The two- loop noise kernel 

For the two- loop noise kernel, the expression analogous to Eq. (|4.12p . from Eq. (I3.24p . is 

X Px('*'^i)/°x(q'^2)/50(k + q, w - - W2) , (4.22) 

where is given by Eq. (I4.10p with decay width given by Eq. (j4.1ip , and is given by an equivalent expression with 
in place of t^. For the range of parameters discussed above, the dominant decay channel for in the thermal bath 
is the two-loop channel (j) — s- 4>XX- The explicit expression for this decay term (and for any of the two-loop terms in 
general) is a very complicated one and there is no simple result available for it, unfortunately. An explicit expression 
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in terms of momenta integrals was derived in [13|, for the on-shell case uu = ti-'(p) only. However, we note that for the 
dominant decay channel (p — > 4>xXi ^(s'*); ^ind since in Eq. (|4.22p we already have a g"^ factor and in each of 

the Px terms we have an extra factor of (from Eq. (j4.11|) ). then A/2 ^ 0{g^h^) and it is subleading compared to 
the one-loop kernel A/i. Because of this, there is no need to proceed further with an analysis of the two-loop term 
here. 

V. CONCLUSIONS 

The effective (microscopic) equations of motion of systems coupled to an environment in the form of a heat bath 
typically include nonlocal dissipation and noise (stochastic) kernels. These terms are generally complicated and need 
to be dealt with somehow for practical purposes for studying the dynamical evolution of these systems. The same 
is true for the evolution of quantum fields in a medium. In this paper we have studied and approached the problem 
of localization of the nonlocal kernels in the context of the dynamics of a scalar background field coupled to an 
environment made of other scalar fields. We have paid special attention to a particular model where dissipation 
emerges as a consequence of a two-stage decay mechanism where a scalar field decays into a light one through a 
catalyst heavy field. This model was identified in Refs. [l^, [lB| and can emerge naturally in the context of Grand 
Unified and Supersymmetric theories, with applications in several contexts, like inflationary cosmology and particle 
physics dynamics in general. 

For this model, we have studied in this paper the nonlocal dissipation and noise kernels that emerge in general 
in the microscopic derivation of the effective equations of motion for quantum fields. From this we have derived 
a fluctuation-dissipation theorem and determined the appropriate conditions for when the kernel terms can be well 
approximated by a local form. 

We should point out that small differences between the local and nonlocal kernels over large evolution timescales 
can still manifest themselves as a secular drift. This, if the scalar field is oscillating, can lead to a phase shift that 
might accumulate, but does not affect the important physical quantities like the rate of dissipation of energy. If the 
dissipation terms dominate over the other kinetic terms, then the overdamped solution emerges and a small correction 
to the equations is equivalent to a shift in the origin of time, which has no physical significance. 

For a scalar field with a expectation value ipc we have found that, in certain situations when the local form is 
possible, the effective equation of motion reduces to 

[□ + M2] ^,(x) + 2r(^, T)^l{x)^,{x) + ^ ^e(x)e(x), (5.1) 

where the stochastic noise source has a correlation function 

(C(x)e(a;')) = 2r{^,T)T5{x - x'). (5.2) 

The situations which we have examined have been for approximately homogeneous fields coupled to the heat bath 
through a catalyst field with coupling g and weak self-coupling A = 0{g^): 

1 . When the catalyst mass is large compared to the temperature and the scalar field is slowly varying on a timescale 
CT~^, where C is a constant which was determined numerically. 

2. When the catalyst mass is small compared to the temperature and the scalar field is slowly varying on the 
relaxation timescale of the catalyst field. 

In case 1, we have also given numerical evidence that the dissipative term is localized on a specified timescale but 
beyond leading order it possesses no analytic expansion with respect to the time coordinate, since beyond this order 
it is dominated by logarithmic terms. 

Finally, we would like to point out a similar model to the one we analyzed here that was studied by the authors 
in Ref. [31|. In that case, they have studied the process of how the homogeneous field (j) can produce fluctuations of 
X through the process of parametric resonance. In that case, oscillations of the scalar field (p lead to the growth of 
the effective mass of Xi which can then cause it to eventually decay into excitations of cr, even if rria is much greater 
than the vacuum mass of cf) or x- Though they analyze a different dynamical regime for the scalar field we think 
that this process of energy transfer from the field (f> to cr, through parametric resonance, is a nice counterpart to the 
quantum process we have studied here. 
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